Simulating accelerated atoms coupled to a quantum field 
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We show an analogy between static quantum emitters coupled to a single mode of a quantum field 
and accelerated Unruh-DeWitt detectors. We envision a way to simulate a variety of relativistic 
quantum field settings beyond the reach of current computational power, such as high number of 
qubits coupled to a quantum field following arbitrary non-inertial trajectories. Our scheme may be 
implemented with trapped ions and circuit QED set-ups. 



Introduction.- The study of accelerated atoms inter- 
acting with a quantum field is a fundamental problem, 
which has attracted a great deal of attention in General 
Relativity. Even though it has been treated extensively 
in the literature [Ull], it still poses intriguing questions, 
as well as experimental challenges for the detection of 
quantum effects induced by acceleration. Moreover, the 
emerging field of relativistic quantum information has re- 
cently increased the attention drawn to the topic, and in 
particular the study of correlations and entanglement in 
non-inertial scenarios j^]- 

Acceleration effects in quantum field theory are well 
understood in the perturbative regime, where, for ex- 
ample, transition rates are obtained for the excitation 
probability of atoms due to general relativity effects. A 
physical paradigm in this regime is the detection of the 
Unruh effect [3| , which, roughly, implies that an atom ac- 
celerated in the vacuum of the field is excited in the same 
way as an inertial atom in an effective thermal field state. 
Going beyond the perturbative regime poses a tough the- 
oretical problem, which gets even harder if one considers 
a set of emitters, in which case we face a many-body sit- 
uation, situation particularly relevant from the point of 
view of relativistic quantum information. 

In this letter we show that a system of accelerated 
atoms [1] coupled to a bosonic field in the discrete mode 
approximation @,[3|, shares interesting analogies to time- 
dependent problems in quantum optics. Our work is mo- 
tivated by the recent experimental progress that has al- 
lowed physicists to develop tools to control the dynam- 
ics of single emitters coupled to fields in set-ups such as 
trapped ions Q or circuit QED |;.9| . The application of 
those systems as analog-simulators of accelerated atoms 
is particularly relevant to many-body non-perturbative 
regimes, where numerical calculations are difficult. Fur- 
thermore, the insight gained by such analogies motivates 
the study of physical effects that may be relevant to the 
experimental detection of the Unruh effect. For example, 
collective phenomena such as superradiance fld\ (which 
are known to amplify the effective coupling of a set of 
emitters to the field) could be used to increase detection 
sensitivity of quantum effects in non-inertial scenarios. 

The letter is structured as follows: We start by present- 



ing the Unruh-DeWitt detector model to characterize an 
atom coupled to a quantum field. Working in the inter- 
action picture from the comoving atom reference frame, 
we show how to account for a uniform acceleration in 
the Hamiltonian. We then proceed to discuss the possi- 
ble physical implementations in trapped ions and circuit 
QED. For both cases a model of emitters with controlled 
time-dependent atom-field couplings is presented, which, 
after some approximations, yields an identical Hamilto- 
nian. For the sake of simplicity we will refer generically 
to single quantum emitters as "atoms" , and the bosonic 
mode as "field" , being the latter a phonon mode in a 
Coulomb crystal of trapped ions, or a photon mode con- 
fined in cavity in circuit QED. Then, as an illustrative 
example, we analyze the physics in the case of a single 
atom to predict the outcome of simple experimental real- 
izations of our ideas. An analogy to a decoupling process 
is explained in terms of the well-known Landau- Zener 
formula. Some future research lines on the topic will 
naturally emerge from our discussion. 

Accelerated Unruh DeWitt detectors.- The Unruh- 
DeWitt detector [J |3, [lH is a standard model for a two- 
level atom coupled to a scalar field. This kind of detectors 
has been extensively used for multiple purposes such as 
acknowledging acceleration effects in cavities, measuring 
quantum correlations between spatially separated regions 
of spacetime, detecting entanglement degradation due to 
the Unruh and Hawking effects and proposing set-ups to 
directly detect such effects. In general, computing time 
evolution under such a Hamiltonian is a complex prob- 
lem, and thus, the affordable calculations reduce to very 
simple scenarios where perturbation theory to the first 
or second order is adequate. 

On the other hand, special interest have the cases 
where the detectors couple only to a discrete number of 
modes. Presented in and extended in ^ such settings 
directly model the interesting case of accelerated atoms 
going through a cavity. Moreover, they also feature a way 
to directly measure the Unruh effect considering cavities 
which are leaky to a finite number of modes 

Let us start by presenting the discretized Unruh- 
DeWitt interaction Hamiltonian consisting in a set of 
two-level atoms coupled to a scalar field. We assume 
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a set of atoms that are accelerating with proper accelera- 
tion a. From the atoms perspective the Hamiltonian can 
be rewritten as 



jm 

(a1^g»["™t(r,«)-fc™2;(r,4)] ^^g-iK„t(r,^)-fc„x(T,^)]-j 

where (r, f ) are the proper space-time coordinates of the 
accelerated detectors, and {t, x) are Minkowskian coordi- 
nates. The following relation holds, 



ct — ^ sinh (ar /c) , x — cosh [ar / c) 



(2) 



Directly from © we see that for constant ^ these co- 
ordinates describe hyperbolic trajectories in space-time 
whose asymptote is the light cone. A uniformly acceler- 
ated observer whose proper coordinates are ([2]), follows 



the trajectory of constant Rindler position ^ = c /a [12 1 



If all the detectors follow such a trajectory we can 
rewrite the Hamiltonian in a form which is suitable to 
find quantum optical analogs. 

Hi = Y,g,„,{a+e^''^^+ H.c.) (aJ„e'*'"M+ H.c), (3) 

jm 

where t is the atoms proper time, and 



-ar I c 



-ar/ c 



a 



(4) 



where we have used the relation Wm = cfcm, and defined 
the parameter a = a/c. We can readily see how we re- 
cover the non-relativistic limit by taking the limit of small 



t pi-f^C^) 
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acceleration, al^e 
Note that the time- independent phase factor, —k„iC^/a^ 
may be simply absorbed in the definition of the particle 
operators. Eq. ([3]) defines the Hamiltonian of interest 
that we aim to simulate. It also reproduces a scenario 
where an array of detectors are resisting near the event 
horizon of a Schwarzschild black hole, if we consider in 
([3]) that a « k/-\/7o, where k is the surface gravity of 
the black hole and /o is the gravitational redshift factor, 
following the arguments in 13 1. 

Physical implementations. Trapped ions.- This setup 
is ideally suited to prepare and measure quantum states 
by well established experimental techniques Q, which 
have found an important ap plic ation in quantum simu- 
lation of many-body physics llj|, and single particle dy- 
namics in special relativity [15|. Previous proposals also 
allowed the simulation of quantum fields using trapped 
ions in general relativistic settings but in different con- 
texts [16|, ll7|, however, paying the price of requiring 
control on the frequency of the qubit and the coupling 
strength. In our proposal, only phase control is needed. 
Our setting is therefore readily exportable to other ex- 
perimental set-ups, as shown below. 

In our scheme we consider a chain of TV ions of mass 
M . For simplicity, we focus on the single- mode version of 



([3]) , and consider that the bosonic mode is a phonon mode 
of the chain, namely, the center-of-mass mode which ac- 
counts to a homogeneous displacement of all the ions. 
State-of-the-art techniques may be used to implement 
phonon sidebands, and to control the atom-field coupling. 

The vibrations of a Coulomb chain consists of a set 
of normal modes described by Hq = ^j^w„ajja„, with 
ijJn the normal mode frequencies, and a„, ajj, phonon 
operators. Levels |0) and |1) are two electronic states 
of the ions. Two sets of lasers couple those levels by 
means of Raman transitions, with amplitudes J^lj , and 
two frequencies WLi.a, (see [8| for details), 

i7L(r) = E ^^M^+e^'^L^.e-^'^^-^-**-*^) + H.c. (5) 

This equation represents a standard atom-light inter- 
action term, with the only peculiarity that we con- 
sider time dependent phases 4>i{t), 4>2{t). Xj are op- 
erators corresponding to the ion displacements, relative 
to the equilibrium positions along the chain. We ex- 
press those displacements in terms of phonon operators, 
Xj = -^j-nXnian+ali), whcre Aij^n are phonon wave- 
functions, and Xn = l/\/2MuJn. 

We choose laser frequencies close to resonance with the 
center-of-mass mode, n = 0, such that ~ — — ^, 
LUh.2 — LJo ~ with f2 <C wq. The coupling ([5|) simulates 
the quantum dynamics of an accelerated Unruh-DeWitt 
detector if: (i) kxn ^ 1 (Lamb-Dicke regime), such that 
we can expand the exponential in powers of Sxj, and 
restrict only to linear atom- field couplings, (ii) ^ 
wq, and /2)kxn <C uiq, |wo — ^n^o\, so that we can 
neglect, in a rotating wave approximation (r.w.a.), all 
couplings to vibrational modes n ^ 0. This yields 



where we have used that the center-of-mass vibra- 
tional modes fulfiUs Mj.o ~ 1/\/N, such that gj ~ 
iVli^.kxol By choosing phases such that ^i(t) — 
~4'2{t) — *I'(t), we arrive at the single- mode version of 
Eq. dS]). Including more modes would require additional 
lasers on resonance with other vibrational modes. Posi- 
tion depending couplings, gj, and frequencies, Q.j, may 
be achieved given a certain intensity profile of the lasers, 
and by focusing multiple lasers on each ion, respectively. 

Typical experimental values are wq, — wqI ~ 1 MHz, 
and kxQ ~ 0.2, such that our requirements are fulfilled 
with values Vti,/2 = 100 kHz, and O = 100 kHz, yielding 
g = 20/-\//V kHz. Those values are well above typical de- 
coherence rates in trapped ions. Also, in order to observe 
analogs to acceleration, <I>(t) has to increase exponen- 
tially over a time-scale comparable to the inverse energies 
involved in the set-up. For example, values a = lO^'^il = 
0.1 kHz, would require to vary the optical phase on times 
scales of 10 ms. This is technically feasible by using, for 
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example, acousto-optical modulators, and standard ex- 
perimental techniques from trapped ion quantum com- 
putation, which indeed require manipulation on a much 
shorter time-scale 18|. Our ideas may also be used with 
spin many-boson models as in [l9| . 

Physical implementations. Circuit QED.- Our sec- 
ond proposed implementation consists of a superconduct- 
ing qubit coupled to a microwave cavity in the strong- 
coupling regime 0. We present a scheme that relies on 
the driving of the qubit frequency only, and is particu- 
larly well suited for our work [20| . The qubits and field 
noninteracting Hamiltonian is 



Ho{r) 



ojQa' a + — 



(6) 



where w is the resonant frequency of the mode, and e 
is the qubit energy. II^{t) describes a driving field, for 
which we assume the following form, 

Hd{T) = - ^ ^ rjj u>d, ,^ cos{uJd^^t + (f)^{T))a2j . (7) 

j u=l,2 

Note that Eq. (O is written as a periodic driving with 
a phase, '/>iy(T), which can be considered to evolve slowly 
in time, in a sense to be quantified below. The qubit- 
cavity coupling in the Schrodinger picture is given by 
Hi — go (f7+ -I- a~) (a -I- a"^). We write the coupling in 
the interaction picture with respect to Hq[t), 

Hi{r) = go5](a+e"^g,(T) +H.c.)(ae-^"«^ -KH.c), 



:,-2i J2v Vsin{uJd^ A'r+'l>v{r)) 



(8) 



where we made the approximation that (I){t) / (j){T) <C 
tOd^-. Consider now that 77 ^ 1, and the choice of fre- 
quencies LUdi j = e — — Wo, J = e — rij + Wo , and 
conditions e,wo, |e — a;o| 3> 50- We expand Q{t) to first 
order in rjj, and keep only slow-rotating terms in a r.w.a: 

Hi{t) « TOo^a+(e'"^-^-*'^i(^)a-he*"^^-*'^^(^)at)-KH.c.. 

3 

This expression takes the form ^ by choosing 0i(r) = 
*(t), 02 (t) = -$(t), and g^ = 50??^- 

In circuit QED, the high energy scales e, ujq are in the 
GHz regime. Low energy scales, such as gj, $7^, may be 
then in the MHz range. Finally, note that <I){t)/4){t) = a. 
Thus, in order to observe effects from effective acceler- 
ation, our scheme requires a on the frequency scale of 
gj, flj, such that the driving fields in JT]) have to be 
controlled with an inverse time in the MHz range. This 
seems technically feasible, since this rate is slower than 
typical evolution times in circuit-QED systems. By using 
different photonic modes in a cavity, as well as local con- 
trol of qubit couplings, the full multi-mode Hamiltonian 
(jSj may be implemented. 

Single detector case: non-adiabatic effects induced by 
acceleration.- We will study a simple case with a two- 
folded purpose in mind: on the one hand it will serve 



to gain some insight on the parameters required for our 
experimental proposal while, on the other hand, it will 
unveil an analogy between non-equilibrium physics and 
quantum effects induced by acceleration. Let us con- 
sider a single atom A, with natural frequency O and uni- 
form acceleration frequency a = a/c = fa^, coupled to 
a single-mode field with proper frequency cj > 17. 

Let us have A excited at t = with no excitations in 
the field. Representing the free eigenstates of atom (A) 
and field (F) with the notation = \AF), the initial 
state would then be \i) = |eO). We let the system evolve 
naturally and concentrate on the probability for A to 
decay to its ground state: 



(9) 



It is convenient to analyze the phase in the Hamilto- 
nian ([3]). The latter is equivalent to a Hamiltonian in 
the interaction picture with respect to a time dependent 
effective frequency 



(10) 



which decreases asymptotically. We can therefore picture 
this situation as a system starting in a certain low energy 
state, in this case matching |e 0) , coupled to a high energy 
state which initially happens to be \gl). The difference 
between the two effective energies varies with time as 



A£;eff(T) = iooe-'^^ - n, 



(11) 



and there will be a level crossing, corresponding to a 
stationary "I>(r), taking place at Tc — \n{ujo / il) / a After 
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FIG. 1. (Color online) Time evolution of VAg for different 
accelerations, with the parameter values g = O.Olfi and ujq = 
1.33f2. The crossing time and the transition probability fit 
the Landau-Zener prediction described in (|12[) . 

a long time (as compared to 1/a) we can obviously ap- 
proximate the low energy level of the system as \g 1), and 
the high as |eO) We can therefore establish an analogy 



4 



between this kind of phenomenon and a typical Landau- 
Zener transition, making therefore a Unk between non- 
equiUbrium physics and quantum field theory in curved 
space-times. If the energy difference happens to vary 
very slowly (which will be the case if a ^ g), we would 
expect no transition between eigenstates to take place, 
so the system will stay in the ground state, which means 
actually decaying into \gl}. If however, the evolution 
happens to be diabatic, the probability of transiting into 
a high excited state (so to say, staying in |eO)) can be 
approximated by the Landau-Zener formula: 



-27rr 



Va 



S( + oo) 



1 



-27rr 



(12) 



where T = / {drAE\r=rJ- Va, is plotted in Fig. □ 

We note several differences with the original Landau- 
Zener model [2l[ . Namely, is not linear but exponen- 
tial in time; also, in our case, only for small couplings we 
have a two-level problem (a r.w.a. approximation cannot 
be performed). Nevertheless, in the regimes considered, 
L-Z is a very good approximation as shown in Fig. [2] The 
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FIG. 2. (Color online) Comparison between Landau-Zener 
theoretical prediction (red dashed) and probability for an ac- 
celerated detector of remaining in the excited state as a func- 
tion of the acceleration (blue solid) with parameter values 
g = O.OISI and luq = 1.33ri. 

atom stays excited until the phase variation rate slows 
down at times close to Tc- Around that moment the cou- 
pling terms get resonant and the decay probability grows 
notably. From then on atom and field get progressively 
decoupled again with the probability stabilizing itself. 

Conclusions.- We have presented a scheme for sim- 
ulating a set of accelerated atoms coupled to a single- 
mode field. First, we have identified the Hamiltonian 
which yields the evolution as seen by a comoving observer 
with the atoms. A method has been presented to ob- 
tain the same Hamiltonian for circuit QED and trapped 
ions by inducing time-dependent sidebands in atom-field 
couplings. Our idea may be extended to many-particle 
experiments which could simulate results that are not 
affordable for classical computers such as arbitrary non- 
inertial trajectories of detectors and many detectors cou- 
pled to quantum fields. Finally, by interpreting our re- 
sults from a Landau-Zener perspective, we have made a 



connection between non-equilibrium Physics and quan- 
tum effects due to acceleration. This reveals an inter- 
esting point that will undoubtedly be exploited in the 
future: general relativistic quantum effects can be stud- 
ied with tools coming from non-equilibrium Physics. 
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